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TITLE 

TYPE 

NUMBER OF WORDS 

TEMPORARY STORAGE 

DURATION 

ACCURACY 

ENTRY 



PRESET PARAMETERS 



Chi-Squared 

Closed with two program parameters 

23 

0, 1, and 2 

(2 Am +3) milliseconds 

+ 2 

When this routine is located at q, entry is made by 



p 


— mF 
50 pF 


p + 1 


26 qF 
00 nF 


p + 2 


any 



DESCRIPTION 



At the end of this routine control is transferred to 
the left hand order at location p + 2. 
When this routine is read in the following preset 
parameters must have been stored in locations 3, k 
and 5, respectively. 

S3 00 F s is a scaling factor and is normally 

00 sF chosen to he 1, 10, 100, or 1000 depending 
on whether none, one, two, or three 
decimal places of accuracy are desired 
after the decimal point. 
S^ 00 F a is the location at which the quantities 
00 aF p., i = 0, 1, ..., m - 1, are stored prior 
to entering this routine. 
S5 00 F b is the location at which the quantities 
00 hF f., i = 0, 1, ..., m - 1, are stored prior 
to entering this routine. 
In its most frequent application chi -squared is given 
by the formula 



~ 2 - 

p m-1 (E - 9. ) 2 
1=0 i 



where each E. is the expected number of members in 
the i " of m classes for a given sample size and the 9 
are the observed values . If we let 



E. 9. 

p. = — S f . = — 2a, b) 

■^l n 9 l n ' ' 

and multiply both sides of equation l) by a number s, 
the resulting equation is 

-*2 m " 1 ^ p i * f lf 
s /. = sn Z . 3) 

This last equation corresponds to the quantity computed 

by this routine. sjL is computed as an integer and 

is placed in the A register at the end of the routine. 

The quantities p. and- f . are fractions and must be in 

the ranges 

< p. < 1 

k a , b) 

< f . < 1 

— l 

A value of 1 for one of the f . may be represented in the 
machine as -1. Each value p. is stored at location 
a + i before this routine is entered, the first address 
a being given by preset parameter Sk, Similarly the 
f . are stored at b + i, b being specified by preset 
parameter S5- The number of values m of the p. or f 
is specified by a program parameter (See Entry). 

The number s is a positive integer specified by 

preset parameter S3 and serves as a scaling factor 

, 2 
for X • Normally s will be chosen to be 1, 10, 100 



- 3 - 



or 1000 depending upon the number of decimal places 

_v 2 

of accuracy required in the value of JL . For example 

ifyt is 2.531* ••• and s is 100, the number in the A 

-39 * 
register at the end of this routine will be 253 x 2 . 

The program parameter n is also a positive integer 
and is chosen so as to put the p. and f. in the required 
range given by 4a, b) as determined by equations 2a, b). 
A logical choice for n is the sample size. The p. are 
then the predicted probabilities and the f . are the 
corresponding observed values The requirements 4a, b) 
will then be satisfied automatically. If the values 
E. and ©. are known directly it may be more convenient 
to chose n to be the smallest power of 10 which is greater 
than or equal to the sample size. In some applications 
a power of 2 may be more convenient. In any case the 
values p. and f . are determined by equations 2a, b), and 
n must be specified upon entering the routine (See Entry). 

If it is much more convenient to produce the values 
p. and f . one at a time, it is possible to enter this 
routine m times with the new values of p. and f . in 
locations a and b, respectively, each time. If this 
procedure is used the program parameter m must always 
be 1, and the necessary summation must be carried on 
outside this routine. Other things being equal, this 
method will be slower and less accurate. 

In addition to 4a, b) there are certain other 

moderate requirements on the quantities involved. The 

first of these is that the following inequality must 

holds p 

y2 m-1 ( p - f f 

JL — = z -J= 1 < 256 = 2 5 5) 

n ■ . p y ' l 

i=0 *i 



* If Library Routine P 1 is used the results may be 
printed in such a way that the position of the 
decimal point is specified. 



- k - 

Since each term in this summation is non -negative, 
each term must also he less than 256. If the latter 
does not hold a division hang -up will occur. If the 
summation is too large the answer -will he in error hy 
a negative integral multiple of 512s. Any danger 
that the requirement 5) might be violated can usually 
he avoided hy using a larger value of n and making the 
appropriate changes in the p. and f . as determined hy 
2a, b). Secondly, in order to obtain the stated 

accuracy of + 2 the product mns should be small 

50 9 
compared to 2 ~ 10 , Although this generally means 

that values of X accurate to a large number of 

decimal places are obtainable, printing out results to 

greater accuracy than actually required or justified by 

the data should be assiduously avoided. 
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1 


42 5L 




Plant address of n 




46 F 




Store m 


2 


L5 19L 
4o 8L 




Set addresses 


3 


46 11L 








l4 f 




Set test constant 


4 


46 20L 








F5 5L 




Plant "link. 
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42 18L 








50 F 




Extract and store n 


6 


00 1043F 

01 20F 






7 


40 F 








41 IF 




Clear Z 


8 


L5 S4 
L0S5 


from 15 


p. - f. 


9 


40 2F 
50 2F 




V 1 1 x 2-8 

P. 


10 


75 2F 
10 8F 






XI 


66 S4 








S5 F' 




Summation 


12 


L4 U 
40 W 






13 


L5 8L 
F4 21L 




Step addresses 
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40 8L 
46 11L 






15 


LO 20L 




Test for end 


16 


36 8L 
50 F 




-39 
ns x 2 in Q 
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75 22L 
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17 


L5 6l 
7^ IF 




s^f rounded off 


18 


00 8F 








26 F 




Exit 
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L5 S4 
LO S5 






20 


75 F 
LO F 
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00 IF 
00 F 






22 


00 F 
00 S3 
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